Ground-state geometric quantum computing in superconducting systems 
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We present a theoretical proposal for the implementation of geometric quantum computing based 
on a Hamiltonian which has a doubly degenerate ground state. Thus the system which is steered 
adiabatically, remains in the ground-state. The proposed physical implementation relies on a su- 
perconducting circuit composed of three SQUIDs and two superconducting islands with the charge 
states encoding the logical states. We obtain a universal set of single-qubit gates and implement a 
non-trivial two-qubit gate exploiting the mutual inductance between two neighboring circuits, allow- 
ing us to realize a fully geometric ground-state quantum computing. The introduced paradigm for 
the implementation of geometric quantum computing is expected to be robust against environmental 
effects. 

PACS numbers: 03.67.Lx, 03.65.Vf, 85.25.Dq 



Oh 
i 

> ■ 
G ■ 

a ■ 

g : 

cr 



> 

in 

(N 
O 
(N 

G\ 
O 
O 



I. INTRODUCTION 

One of the major issues in the present quantum infor- 
mation technology is the protection of quantum systems 
from external noise. The noise can be due to impre- 
cise control or undesired interaction with the environ- 
ment. To eliminate some errors of the first kind, geo- 
metric quantum computing (GQC) has been proposed 
[ll-Q. The idea behind it is to use the so-called geomet- 
ric phases Q to manipulate the quantum state. Usually 
the appearance of a geometric phase is associated with 
a cyclic and adiabatic manipulation of the system. The 
advantage is twofold: the adiabatic evolution eliminates 
the problem of fine-tuning of the evolution time and the 
geometric operators can be intrinsically robust against 
fluctuations of the imprecise control field @. 

It is not possible to realize universal quantum comput- 
ing only by means of the Abelian Berry phase Q . There- 
fore, utilizing the non- Abelian geometric phases Q, the 
geometric quantum computing was proposed in 1999 [l|. 
In this case, the logical states are encoded in a degen- 
erate subspace of the full Hilbcrt space and they can be 
manipulated by means of purely geometric non- Abelian 
operators. The first proposals of the physical implemen- 
tations were based on trapped ions and exploited the so- 
called tripod Hamiltonian to build a doubly degenerate 
logical subspace [H, . Subsequently, the same Hamilto- 
nian has been used in many quantum systems providing 
means for possible implementation for neutral atoms, su- 
perconducting circuits, quantum dots, and Bose-Einstein 
condensates |10l - tl4j . In conjunction with these proposals, 
various studies have been carried out to test the stability 
of geometric gates against different types of noise. Here, 
GQC has proven robust against fast adiabatic fluctua- 
tions of control fields [l5| and certain type of environ- 
mental noise fl6l . [I?} ■ 

Despite the theoretical efforts, we are still far from a 
full GQC, the experimental observations of non- Abelian 
phases remain very limited and still under discussion [l8| . 



One possible reason behind the absence of experimental 
verification is the decoherence due to the interaction of 
the quantum system with the environment. In fact, be- 
cause of the condition of adiabatic evolution, it is chal- 
lenging to protect the system for long enough time to 
implement geometric operations. This is the case for all 
of the above proposals based on tripod-like Hamiltoni- 
ans in which the evolution does not occur in the lowest- 
energy eigenspace. In general, one expects that decoher- 
ence such as relaxation can be avoided if the evolution 
occurs in the ground-state manifold. The possibility to 
exploit this robustness in large systems to perform adia- 
batic quantum computation is still under discussion fl9| . 
However, the robustness of ground state evolution has 
been recently confirmed for a driven two-level system in 
the adiabatic limit [2(il - l22j . In the same studies, it was 
shown that relaxation induced by the environment can 
help to keep the system in the ground state potentially 
leading to an improved robustness of the ground-state 
evolution. 

These results suggest that GQC which exploits a de- 
generate ground-state eigenspace can be experimentally 
feasible. The first step in this direction has been taken in 
Ref. [23[ in which a way to detect non- Abelian phases in 
superconducting Josephson devices was proposed. How- 
ever, the work mainly focused on observing the geometric 
non- Abelian effects and only a particular unitary trans- 
formation was discussed. 

The extension to a universal set of transformations is 
non-trivial and is the main result of this paper. Starting 
from a general Hamiltonian, we show how to obtain the 
basic logical gates in an abstract context. Then, using 
the same system as in Ref. [23| , we show how these gates 
can be implemented in a physical set-up. We build all the 
logical gates needed for quantum computing: a complete 
set of single-qubit gates and a two-qubit gate exploiting 
the mutual inductance between two circuits next to each 
other. Together these gates form an universal set allow- 
ing full geometric quantum computing. 
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The paper is organized as follows. In Sec. [H] we in- 
troduce a model Hamiltonian with a doubly degenerate 
ground state which can produce universal single-qubit 
transformations. In Sec. IIII1 we present the physical 
system and show how to map the abstract model into 
the physical Hamiltonian. In Sec. IIV| we show how to 
manipulate the quantum system in order to obtain two 
fundamental single-qubit gates. Section [V] is devoted to 
the implementation of a geometric two-qubit gate. In 
Sec. I VII the readout measure is discussed and Sec. IVIII 
concludes the paper. 

II. MODEL 

We consider a quantum system with a two-fold de- 
generate ground-state eigenspace and an excited state. 
The system is steered adiabatically along a loop by ex- 
ternal control fields. If the evolution is adiabatic and 
occurring in time T a d, the dynamics can be accurately 
determined by instantaneous diagonalization of the sys- 
tem Hamiltonian since the system follows approximately 



the evolution of the eigenstates. It is convenient to fix 
the logical basis in order to have the initial Hamiltonian 
in a diagonal form: in particular, the two degenerate 
ground-states used as logical states are denoted as |0) 
and |1), and the excited one representing an auxiliary 
state is denoted as \a). Taking the ground and excited 
state energies equal to 77 and 77 + a, respectively, we have 
H(t = 0) = t]l+a\a)(a\, where 1 = |0)(0| + |l)(l| + |a)(o|. 
The basis {\a), |0), |1)} is time- independent and corre- 
sponds to the eigenbasis only if the Hamiltonian is in the 
above form. 

We require that the Hamiltonian has always a doubly 
degenerate ground state during the control cycle. Thus 
we can formally write it as H{i) = r/1 + a|t>i)(ui| where 
|t>i) is a linear combination of the time- independent 
{\a), |0), |1)} states. Throughout the paper we adopt 
the convention to denote with the operators with hat 
symbol (e.g., H) and their representation in the basis 
{|o>, |0), |1)} with the tilde (e.g., H). Following Ref. [H, 
the most general representation with the spectrum of H 
in the basis {|a), |0), |1)} reads 



H 



i] + a cos 2 9\ cos 2 62 —ae cos0i cos 2 82 sin0i —ae i( ^ 2 cos0i cos #2 sin #2 

-ae* cos 6\ cos 2 82 sin 8\ i] + a cos 2 82 sin 2 8\ at 
-ae* cos 0i cos 2 sin 2 ae - ** - * 2 ) cos0 2 sin0i sin0 2 rj 



- a sin' 



OS V\ COS t>2 sin t>2 

cos 82 sin 8\ sin 82 



(1) 



where 0i, 02, 81, and 82 are experimentally modulated 
parameters. In the following we will assume that 0i = 0. 
The parameters 8\, 82, and 02 describe effectively the 
time evolution of the Hamiltonian. The instantaneous 
eigenvectors of H can be written as 

\vi) = -e~*cos0icos0 2 |a) + e~* cos0 2 sin0i|O) 

+ sin 8 2 1 1), 
\v 2 ) = sin 6»i I a) + cos0i|O), 

\v 3 ) = e"* cos 6»i sin 6» 2 1 a) - e"^ 2 sin0i sin0 2 |O) 



-cos 6» 2 II) 



(2) 



where ^2) and l^) are the degenerate time-dependent 
ground states. 

The restriction to adiabatic evolution, i.e., aT ac i/h S> 
1, assures that no transitions to the excited state occurs; 
the evolution is restricted to the ground-state eigenspace. 
The final unitary transformation depends only on the ge- 
ometric features, i.e., on the loop 7 covered in the param- 
eter space {6*i, 8 2 , 4> 2 }, and it reads 0, 0, [H| 



(A ei d8i + Ae 2 dd 2 + A^dfo) , (3) 



where V denotes the path ordering operator. The A x 
are the so-called connections and, restricted to the de- 
generate eigenspace, they are defined by A^—i i—i = 



(vj\d x \vi) = 2,3). In our case, restricted to the de- 
generate eigenspace, they assume a simple form 



A 6l = 
Ae 2 = 

A(h = 







sin o 2 



3 -i</>2 



sin6» 2 









-i sin 2 02 



(4) 



Since Ag 2 gives no contribution, the unitary transfor- 
mation associated with a loop 7 in the parameter space 
{0i, 02, <fo } can be written as 



Vcxp 



(A 6l d9 1 +A 4 



(5) 



III. PHYSICAL REALIZATION 

The physical set-up considered for the implementation 
of GQC is shown in Fig. QJa). It is composed of three 
superconducting quantum interference devices (SQUIDs) 
in series with two superconducting islands between them 
[23|. The SQUIDs are operated as tunable Josephson 
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FIG. 1: (Color online) (a) Scheme for the superconducting geometric qubit. The gate voltage and the capacitance of the ith 
island are denoted by V 9i and C gi , respectively. The ith SQUID is described by Josephson energy Ei, capacitance d, and 
flux The phase across the device is (p. (b) The hexagonal (honeycomb) stability diagram in the n gs — n g d plane with all 
the SQUIDs closed (Ei = 0). In the inset, the loops that generate the logical gates are provided. The trapezoidal loop (solid 
line 1— > 2 — > 3 — > 4 — > 1) generates the Hadamard gate and triangular loop (dashed line 1 — > 4 — > 5 — > 1) generates the U2 
tranformation. The circle represents the point in which the SQUIDs are closed and the phase across the device is reversed from 
p to — ip in the implementation of the U2 tranformation. Notice that away from the degeneracy lines the SQUIDs are open 



(Ei 7^ 0) and the loops represent the projection on the 



rigd plane of the degeneracy curves in the whole {n ga ,n g d, Ei} 



space. Here, we have chosen k = 1/3 and a big energy gap a = Ec to visualize better the loop. 



junctions which can be closed (Josephson energy Ei = 0) 
and opened (Ei ^ 0) by controlling the magnetic flux 
through them. The phase difference of the order param- 
eter across the whole device, ip = <p R — 4>l, is controlled 
by the magnetic flux <J> through the outer-most loop. A 
slightly simpler device has been already used to measure 
experimentally the Berry phase [25| . 

The total Hamiltonian -ff p h ys is composed of the charg- 
ing Hamiltonian, if c h, and the Josephson Hamiltonian, 
Hj. The first one reads 



Hch =E Cl (fii - n gi ) 2 + E C2 (n 2 - n g2 ) 2 
+E m (hi - n gi )(h 2 -n g2 ), 



(0) 



where fii is the operator for the excess number of Cooper 



pairs on the island i and 



C gi V g J(2e) 



is 



the 



corresponding gate charge 
E Cl = 2e 2 C E ,/C 2 , E C2 
ie 2 C m /C 2 . 



The charging energies are 
2e 2 C Ei /C 2 , and E m = 
Here, C^. is the total capacitance of the 



island i, namely, 



G 



Cr + C m + C„„ + 



— Cl + C rn 
Cl, and C 2 



^C gi 



E 



Cf and 

r, 2 



Above, the C° is the self-capacitance of island k. 
The Josephson Hamiltonian is given by 

^ 00 

H.i = -^ J2 (^L($i)e W2 |n 1 + l,7i 2 )(7 11 ,n 2 | 

Tii ,ri2 — — 00 

+E m ($ m )\n 1 + l,n 2 - l){ni,n 2 | 
+E R (<S> R )e-^ 2 \m , n 2 + 1) (m , n 2 | + h.c.) , 



(7) 



where \ni,n 2 ) denotes the eigenstate of h\ and h 2 and 
{Ei(&i)} are the tunable Josephson energies which are 
controlled by the fluxes through the SQUIDs {$i}. 
Above we have neglected the contributions arising from 
the finite loop inductance and suppose that the SQUIDs 
are perfectly symmetric. 

It is convenient to use the total charge on the two is- 
lands, h = hi + h 2 , and the charge asymmetry, rh = 
h\—h 2 as quantum numbers. Thus, the quantum state of 
the system is denoted by \n, m) and the new gate parame- 



ters are n gs 



L gi 



'U-2 



and 



n gd 



. If the charg- 



ing energies of the islands are equal Ec = Ec 1 = Ec 2 , 
the c harg ing Hamiltonian © with the new notation 
reads 



2G] 



H ch = 



2E C + E m 



1- \2 , ^ E C ~ E m 2 

(n-n gs ) +— — —£-{rn-n g d) 



2Ec + E m 



(8) 



and the Josephson Hamiltonian can be expressed as 
1 °° 

Hj =-- ^ (E L (<$>L)e lv/2 \n+\,m + l)(n, 



in 



n.m— — 00 



+ E m (<f> m )\n,m + 2)(n,m\ 

+ E R ($ R )e- lv / 2 \n+ l,m- l)(n,m\ + h.c.V (9) 

If all the SQUIDs are closed, i.e., E L = E m = E R = 0, 
the conventional stability diagram with a hexagonal lat- 
tice structure in the n gs — n g d plane is recovered as shown 
in Fig. [TJb) . Hexagonal cells of the lattice are de- 
noted by h nm and inside the cell the state |n, m) is the 
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non-degenerate ground state. The edges of the cells rep- 
resent lines of double degeneracy and the intersection 
points of three cells are associated with triple degenera- 
cies. 

In the following, we restrict the evolution to the vicin- 
ity of the triple degeneracy point {n gs ,n g d} = {(3 — 
&)/2,0} with k = (2E C - E m )/(2E C + E m ) at which 
the states \n = l,m = 1) = \a), \n = l,m = —1) = |0) 
and \n = 2, m = 0) = |1) have the same charging en- 



ergy They correspond to 1 excess charge on the first 
island and no charge on the second, no charge on the 
first island and 1 in the second, and 1 charge on the 
first and second island, respectively. For small Joseph- 
son coupling, Ei -C Eq, and in the vicinity of the triple 
degeneracy point, the system evolution is restricted to 
the basis {\a), |0), |1)}. Thus we can truncate the Hamil- 
tonian iJ p i iys into this basis with the representation 



/ E c [(l-n gs ) 2 + k(l-n gd ) 2 ] -\E m ($ m ) 



phys 



\ -\E R & R )e-* 



Ec[(l 



i ) 2 + k(l + n gd ) 2 



-i J B L ($ L )e^ 



-\E R {^ R )e^ 
-\E L (<S> L )e- l -% 
E c [(2 - n gs ) 2 + kn 



(10) 



gd\ 



where E c = (2Ec+E m )/i. 

The mapping between the matrices H and -ff p h ys is 
straightforward for the diagonal elements and the ab- 
solute values of the off-diagonal terms. For example, 
we can solve the system of equations for the diagonal 
terms Ha = -ff p hy S ,M- From the first equation we obtain 



V = E C [{1 



,) + fc(l — n gd ) 2 } — a cos 2 8 X cos 2 



we solve the remaining equations to have 




, 3 cos" 



4£c 
cos(26>i) cos 2 6 

1 4kE c 



+ 2(3-fc) 



and 



(11) 



Mapping the absolute values of the off-diagonal terms 
as a function of 8\ and 02 results in 



\E m {$ m )\ 

\E R ($ R )\ 



= 2ol cos 62 sin 8\ sin 82 
= 2a cos 81 cos 2 62 sin 0i 
= 2a cos 81 cos 82 sin 62- 



(12) 



The parameter a is the energy gap of the system but it 
has no general simple expression in terms of the physical 
parameters. In the following, we impose that a is con- 
stant during the evolution and determine the correspond- 
ing evolution of the physical parameters by Eqs. (fTT|) 
and ([T2"|) . However, the value of a and even its time- 
independence are not important for the model since the 
geometric operator does not depend on the value of the 
energy gap. This further degree of freedom can be uti- 
lized to design a timc-dcpcndcnt a leading to ease the 
control of the physical parameters for the experimental 
implementation. 

Even if Eqs. (fTTj) and (TT2|) arc satisfied, comparison 
between Eqs. ([1]) and (|TU| reveals that it is not possible 
to match the phases of the off-diagonal terms, and thus, 
a global bijective mapping does not exist. In fact, if all 
the SQUIDs are open, it turns out that the two Hamil- 
tonians have different spectra and the physical Hamilto- 
nian has a degenerate ground state only for particular 
choices of the phase across the device: <p = 0, ir (see also 



J26|). However, there is still the possibility to map the 
two Hamiltonians locally if we assume that it is possible 
to close the SQUIDs perfectly. 



IV. SINGLE-QUBIT LOGICAL GATES 

Since the calculation of the geometric unitary trans- 
formation is simpler using the model Hamiltonian H, we 
use the following approach: first determine the loop for 
H which produces the desired geometric transformation 
and then map it into the physical parameter space. 

Along a single loop it is convenient to keep (f>2 constant 
thus restricting the curve to the 8\ — 82 plane. It follows 
from Eqs. (U]) and (O that a rectangular loop in the 81 — 82 
plane produces the final transformation (in the logical 
{|0> 5 |1}> basis) 



r 



cxp 







e' 

2 9 



where g = (sinflf - sin^ l )(6»f / - 0f ), and 8f [B, 
the maximum (minimum) value of 8i . 



(13) 



m ) is 



A. Hadamard gate 

We set <f>2 = and perform the loop 

(0i,0a) : (0,0)->-(0,7r/6)^(7r/2,7r/6) 

-►(tt/2, 0)^(0,0), (14) 

which corresponds to S^ 1 = tt/2, 8^ = 7r/6, and 0™ = 0. 

The corresponding transformation using Eq. (|13l) is the 
Hadamard gate (apart from an irrelevant global phase), 



J_/l -1 
V2K 1 1 



(15) 




FIG. 2: (Color online) Time evolution (normalized by the adiabatic time T a( j) of the physical parameters to implement the 
Ui transformation, (a) The gate parameters n gs (solid line) and n g d (dashed line) and (b) the SQUID parameters \E m \ (solid 
line), \El\ (dashed line), and \Er\ (dot-dashed line) in units of Ec- The blue dot denotes the point, in which the SQUIDs 
are closed and the phase across the device is reversed tp — > —tp. In all the plots the energy gap is constant a = 0.1 x Ec and 
k = 1/3. 



In this case, the mapping to the physical Hamiltonian 
is simple. If we set the phase tp across the device to 
zero, the Hamiltonian corresponding to Eq. (|10p can be 
globally mapped into the one in Eq. ([1]) with the help 
of Eqs. (fTTj) and (fT2j) . The corresponding loop in the 
rigs — n g d plane is shown in the inset of Fig. [TJb). 

B. Phase gate 

Let us consider the case in which the model Hamilto- 
nian ([1]) undergoes a loop with time-independent ^ / 0. 
We choose Of = it/2, Of = tt/2 and the loop 

(0i, 2 ) : (0, 0) -> (tt/2, 0) -> (tt/2, tt/2) 

^(0,^/2)^(0,0), (16) 

which produces the transformation 

The choice of this particular loop is due to the sim- 
plicity in calculating the unitary operator. In fact, the 
path (tt/2, tt/2) — > (0,7r/2) is not necessary since all the 
SQUIDs are closed (Ej = 0) and the n gs and n g d are con- 
stant. The Hamiltonian ([1} does not change along this 
path [see Eqs. (fTTj) and (fT2|) ] which therefore is physically 
irrelevant. 

To locally map the physical Hamiltonian H p hys to H 
along the loop in Eq. (|16[) . the physically irrelevant step 
(tt/2, tt/2) — > (0,7r/2) can be neglected. The physical 
evolution can be obtained from Eqs. (fTTj) and (fT2]) and 
it is shown in Figs, [ljb) and [5J The sequential open- 
ing of SQUIDs R, L, and m in Fig. &b) corresponds to 
the three edges of the triangular loop in Fig. [ljb) . How- 
ever, particular attention must be paid to the choice of 
the phase across the device tp for the correct mapping 
between the Hamiltonians. 



In the beginning of the loop all the SQUIDs arc 
closed and we fix the phase to tp = 2<p2- At point 
(#1,02) = (tt/2, 7r/2), again all the SQUIDs are closed, 
both the Hamiltonians are diagonal and the mapping be- 
tween them is trivial. At this point [marked with a dot 
in Figs. QJb) and[2[b)], no information about the phase 
is present, and hence we can adiabatically reverse the 
phase across the device tp to —2<f>2- During this pro- 
cess, an additional diagonal contribution proportional to 
tp arises in the physical Hamiltonian due to the ac Joseph- 
son relation. It appears in the Hamiltonian as a finite 
bias voltage and can be treated as an effective shift in 
the gate charges 8n gi in Eq. and hence, it can be 
eliminated by an proper choice of n gi allowing to restore 
the desired degeneracy. After this procedure, the right 
SQUID is open but, since the phase through the system 
has been changed, it is possible to map the Hamiltoni- 
ans along the third path. The final result is a mapping 
between the physical Hamiltonian and the model Hamil- 
tonian with constant phase <fo- Thus, this generates the 
U2 (4>2) transformation. 

Having established the loop to build a logical U2 gate, 
we can obtain a phase gate if we perform two sequen- 
tial loops with opposite phases 4>2 &ud — 4>2- The final 
transformation is 

c/ ph (0 2 ) = u 2 (-4> 2 ) u 2 (0 2 ) = % 2 ) - 

(18) 

which allows us to control the relative phase between the 
logical states. 

C. Sequential application of logical gates 

To obtain the logical transformation JJ\ and C/ 2 ((fo) 
we have assumed that <p2 is constant. However, in a 
sequential application of different logical gates we must 
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change 4>2 from one loop to the other. These changes have 
no effect on the dynamics of the system. In fact, they 
are performed when the SQUIDs are closed, 0\ = 02 = 
and thus, from Eq. Q even if changes the associated 
connection vanishes. 



V. TWO-QUBIT LOGICAL GATE 

To complete the set of gates for universal geometric 
quantum computing we need a non-trivial two-qubit gate 
[27J. To this end, we use two superconducting qubits as 
described above. To couple the qubits, the mutual in- 
ductance between the two circuits is exploited as the key 
physical phenomenon as depicted in Fig. [3j If a current 
is flowing through one circuit, the flux and the overall 
phase across the other one are changed allowing us to 
implement controlled logical operations. In the following 
we use the first circuit as a control qubit and the sec- 
ond one as a target qubit. Thus, we let the current flow 
through the first circuit and apply to the second qubit 
a geometric transformation which depends on the flux 
induced by mutual inductance. 



A. Current through the first circuit 

To couple the qubits, a stationary current is induced 
through the first circuit whose direction depends on the 
logical state of the qubit. 

Denoting the Hamiltonian associated with the first cir- 
cuit as -ffphys, 1 and the relative phase across it as ipi, the 
current operator throu gh the first circuit is defined as 
I\ = (2e/ h)d Vl H p hys,i UM- To calculate the current at 
ipi = 0, we consider the variation of the Hamiltonian 
for ip\ <C 1: -ffphys, i(0) +&Hph y8 ,i(v?i)- The Hamiltonian 
-ffphys, i(0) is the dominant contribution and 8 -ffphys, 1(^1) 
is a perturbation that can be obtained as linear contribu- 
tion in 931. The dominant part -ffphys, i(0) can be mapped 
into H and analytically diagonalized. Thus, the eigenba- 
sis in Eq. ([2]) represent the preferred basis for our calcula- 
tion. Using the degenerate-state perturbation theory, it 
is possible to diagonalize 5 H p hy s,i(vi) m the degenerate 
subspace. The eigenstates are \v' 2 ) = (i|ua) + \ V 3))/V% 
and \v' 3 ) = (— i\v2) + \v 3 ))/V2 and, in the new basis 
j>i), \v' 2 ), \v' 3 )}, using Eq. {IT} and (JT2J) , Hph y8 ,i(</?i) can 
be written as 



-ffphys.lC^l) 



T] + a 
-ffi V 



-r<pi 

a cos 2 82 sin (28i ) sin 6*2^1 




f<Pi 


rj + a cos 2 82 sin (20 1) sin^v?! 



(19) 



where / = a [i cos (20i) + sin (20 x ) sin6» 2 ] sin (20 2 ) 1(2^2). 

Using Eq. (|19[) . the current operator I\ for (p\ = 
can be calculated and the average current associated 
with the degenerate eigenstates reads = = 
(-l) i+1 (2ea/ft)sin(26»i)sin6i2Cos 2 6'2 (with i = 2,3). 
Thus, the current flows in opposite direction depending 
on the system state. 

However, this is not sufficient. In fact, if initially the 
SQUIDs are closed, the current states \v' 2 ) = l/y/2(i\0) + 
|1)) and \v' 3 ) = 1/a/2(— *|0) + |1)) do not coincide with 
the logical states. Since we want the logical states to be 
associated with the flowing current, we must rotate the 
logical basis into the current state basis before opening 
the SQUIDs using the geometric operators U p h(ir/2)Ui. 
After this rotation, the opening of the SQUIDs could it- 
self produce a mixing of the \v' 2 ) and \v' 3 ) states. To avoid 
this we open the SQUIDs following the path (01,02) ■ 
(0,0) -> (0,7r/4) -> (7r/4,7r/4). Using equation simi- 
lar to Eq. ([5]), by direct calculation, it can be verified 
that the current states \v' 2 ) and \v 3 ) are not mixed along 
this path [28[ . Thus, the SQUIDs are open, the current 
flowing through the first circuit is I\ = ±ea/(hV2) de- 
pending on the initial logical state |0) or |1). 

Equation (flT))) allows us to estimate the effect of self- 
inductance L on the system. Even if the phase across the 



device is initially set to zero, the current I\ flowing in the 
circuit induces a change <5$i in the magnetic flux. Thus, 
the system experiences a small additional phase 5ip\ = 
7n5<f>i/$o = ±7rL|/i|/<f>o associated with the eigenstates 
\v' 2 ) and \v' 3 ). As observed from Eq. (fl9|) . because of 
the different sign of the phase shift, both the degenerate 
eigenstates have the same energy and the degeneracy is 
preserved up to the linear order in Sipi. 



B. Two-qubit interaction Hamiltonian 

Let us assume for simplicity that the phase across the 
second circuit, f 2l is initially fixed to zero. Due to the 
current through the first circuit, the second one expe- 
riences an additional flux through the outermost loop 
<S<I>2 = ±M2i, where M is the mutual inductance be- 
tween the loops. This induces a phase increment 5<f2 = 
118^2/^0 = ±7rA/|7i|/$o depending on the initial logical 
state of the first qubit. 

The total Hamiltonian is -fftot = -ffphys, i(</ 3 i = 0) + 
-ffphys, 2(8^2)- We suppose that the quantum state of the 
first circuit is a combination of the degenerate states \v' 2 ) 
and \v 3 ). In this subspace, the current operator I\ is 
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diagonal and we can write the phase operator 5if>2 as 



Sif> 2 



7TM\h\ 
$ 



(\v' 2 ){v' 2 \-\v' 3 ){v's\)- 



(20) 



In the second circuit we suppose to have always at least 
one SQUID closed and thus, no current flows through it 
and the first circuit does not experience an additional 
mutual flux. The charging Hamiltonian of the second 
circuit does not depend on the phase operator 5(f>2- For 
this reason, the interaction between the two circuits oc- 
curs through the Josephson Hamiltonian of the second 
2 ). Using Eq. © it can be written as 



h 



-h 



circuit H^\{5(p 



-- ^ {E L ^LaV^ ® \n + l,m + l)(n,m\ 

n,rn 

+^m,2($m,2)l ® \n, TU + 2) (n, 77l| 

+E R:2 {<S>R,2)e~ t ^ L ® |n + l,m- l)(n,m| 



M 



+h.c 



[ J,2 



OV2) 



|00) -> e^ 2 |00), 

1 10) -> e-^ 2 |10), 

101) -> e-^ 2 |01), 

|11) -> e l ^ 2 |ll). 



(23) 



(21) 



where -Efc,2 and indicate the Josephson energy and 
the flux of the fcth SQUID of the second circuit, respec- 
tively, and S(p2 acts on the Hilbcrt space associated with 
the first circuit. 

Writing the Hamiltonian in the degenerate subspace 

{K) 5 K)}> we have 



$ 2 ± MIj 



(22) 



where Hj,2(iSip2) is now in the form of Eq. (0). This 
Hamiltonian allows us to perform an operation on the 
second qubit depending on the state of the first one. In 
fact, if the first qubit is in the state \v' 2 ) (\v$)), we act on 
the second by the Hamiltonian Hj^iS^) [Hj^—S^)]- 



C. Two-qubit logical gate 

To be more specific, we choose to implement a par- 
ticular controlled phase gate on the second qubit. The 
SQUIDs are initially closed and the degenerate space is 
spanned by the logical states {|00), |01), |10), |11)}. We 
rotate the states of the first qubit by the transformation 
J7ph(f/2)?7i- In this way, when the SQUIDs are open the 
current flows through circuit depending on the initial log- 
ical state (for example, for |0) we have current 1\ and for 
|1) we have current — The total phase across the sec- 
ond circuit is ±5<f2 and we apply a phase gate discussed 
in Sec. HVBl Then, we close the SQUIDs of the first 
circuit and apply a transformation [f7 p h(7r /2)U\y back 
to the logical basis. At the end of this procedure, the 
logical states acquire a phase proportional to 5<f2 



FIG. 3: (Color online) Circuit scheme for two qubit gate im- 
plementation. The circuits are coupled by mutual inductance, 
M. A current ±7i circulating in the first circuit corresponds 
to an additional magnetic flux S<&2 = ±MJi through the sec- 
ond. 



The total transformation can be written as 

E/ 2q ub s = e***»°-' ® |0)(0| + e -«w ® |1)(1| 
= ( e iS ^ a - ® 1) 

x(l<g> |0)(0| + e - 2 *<W* <g, \1){1\), (24) 

where a z = |0)(0| - |1)(1|. 

From the last line in Eq. (|24|) . we observe that the 
operator represents a controlled phase operation. An ar- 
bitrary controlled phase gate can be produced increasing 
the current I\ in the first circuit or by iterative applica- 
tion of the presented logical gate. The additional single 
qubit <j z rotation can be compensated after the iteration. 



VI. MEASUREMENT SCHEME 

The logic information is encoded in the charge degree 
of freedom of the two islands and the logical states cor- 
respond to |m = 0, n<x = 1) = |0) and \n\ = l,n 2 = 
1) = |1). Thus the measurement of the charge on the left 
island is sufficient to distinguish between the two logical 
states. This measurement can be done in several ways, 
for example, by using a radio frequency single-electron 
transistor (rf-SET) [29l - l3l| that is capacitevely coupled 
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to the left island. The SET is embedded in a tank cir- 
cuit and the readout of the charge state is performed by 
monitoring the reflected radio frequency signal from the 
circuit sensitive to the charge state on the island. 

This measurement scheme has been proved to be fast 
and sensitive |3(|. An additional advantage is that the 
measurement apparatus is switched off when no radio 
frequency signal is induced in the resonant circuit. This 
point is critically important in our case since the inter- 
action between the system and an external measurement 
apparatus can lead to the breaking of the ground state 
degeneracy. The effect of the additional measurement 
circuit on the island increases its total capacitance and it 
can be taken into account by a compensation in the gate 
voltages. 

VII. CONCLUSIONS 

In this paper, we have presented means to implement 
ground-state geometric quantum computing in supercon- 
ducting systems. We have shown how to build all the log- 
ical gates necessary for a geometric quantum computing 
by manipulation of a few experimental parameters. The 
main advantage with respect to the previous proposals 



for geometric quantum computing is that the evolution 
occurs completely in the ground state. This feature ren- 
ders the implementation robust against low-temperature 
environments which motivates for an experimental verifi- 
cation. In fact, the physical system discussed here can be 
potentially built and similar measurements on a simpler 
system have been already carried out p5| . 

Although we focused on the implementation in super- 
conducting circuits, it is possible that similar Hamiltoni- 
ans can be found in other physical systems as happened 
for the first experimental proposal in Refs. [1, 

The future challenge is to study in detail the robust- 
ness of our scheme against environmental effects. In this 
context, the results of Ref. [2(il - [22j seem promising since 
the non-degenerate ground state evolution has proven to 
be robust against low-temperature environmental noise. 
If similar results are obtained for the present model, it 
could become a new paradigm for the implementation of 
robust geometric quantum computing. 

We acknowledge Academy of Finland and Emil Aalto- 
nen Foundation for financial support. We have received 
funding from the European Community's Seventh Frame- 
work Programme under Grant Agreement No. 238345 
(GEOMDISS). 
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